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Abstract. Magnetic fields probably play a central role in the dynamics and 
thermodynamics of ICMs and their interactions with AGNs, despite the fact 
that the fields usually contribute relatively little pressure; i.e., the ICM is a 
"high-/?" plasma. More typically, the roles of magnetic fields come through 
"microscopic" influences on charged particle behaviors, and through magnetic 
tension, which can still be significant in subsonic, high-/3 flows. I briefly review 
these issues, while exploring the underlying question of using the commonly- 
applied magnetohydrodynamics model in the ICM when Coulomb scattering 
mean free paths can sometimes exceed tens of kiloparsecs. 



1. Introduction 

Powerful AGN jets typically deposit most of their momentum and energy in low 
density, fully ionized and largely collisionless plasmas that constitute the intr- 
acluster media, or ICMs. Inevitably in such media, charge mobility differences 
lead to electrical currents and, thus, to magnetic fields. Those magnetic fields 
limit the motions of the charged particles, affecting the momentum, energy and 
charge transport through the plasmas and the electric currents establishing the 
magnetic fields themselves. The detailed physics of these interactions is com- 
plex and most accurately explored on "microscopic" levels through the tools 
of plasma physics. On the other hand the vast degrees of freedom inherent in 
plasma treatments, especially when applied on "macroscopic" scales, often makes 
such treatments unwieldy. It is common, instead to model these media, includ- 
ing their interactions with AGN jets, through magnetohydrodynamics (MHD). 
Like other continuum approximations, MHD carries with it assumptions that 
can potentially mask or exclude relevant physics. On the other hand, when ap- 
propriate, MHD provides a very powerful tool, and sometimes the only practical 
tool for the exploration of the interactions at the center of discussion in this 
meeting. It turns out that even nominally weak magnetic fields can have pro- 
found influences on both flow dynamics and thermodynamics in these settings. 
I was asked to explore these issues in this presentation. 

I begin with a short review of the assumptions built into MHD and an eval- 
uation of their applicability in ICMs. Concluding, with some caveats, that the 
MHD model is appropriate on many important length and time scales in clusters 
I summarize some common dynamical and thermodynamical MHD issues and 
then discuss some specific AGN/cluster interactions where the presence of mag- 
netic fields appear to be important. My brief comments are necessarily incom- 
plete. More thorough discussions of MHD and its connections to plasma physics 
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are widely available (e.g., Bovd Sz Sanderson 20031 : iGoedbloed fc Poedts 20041 : 
Kulsrud 200llPriest Forbes 2000l )~ 



I do not consider the origins of magnetic fields in clusters or radio galaxies, 
although that is of much current interest. The proceedings of the conference, 
"The Origin and Evolution of Cosmic Ma gnetic Fields", pr ovide a good intro- 
duction to many aspects of the problem (|Beck et al. 200(3 ). Suffice it to say 
that both observational and theoretical estimates of magnetic field strengths 
in clusters are typically in the vicinity of a few microGauss , give or take an 
order of magnitude (e.g., Carilli Taylor 20021: Dolag 200^: iKronberg 19961 : 
Ryu et al. 20071 : [Schekochihin et al. 20051 : IVikhhnin et al. 200ll ): similar, if some- 
what larger, estima t es apply in the lobes of FRI and FRII radio galaxies (e.g.. 



Croston et al. 20051 : Worrall &l Birkinshaw 2000l ). I will take a microGauss for 



a characteristic magnetic field strength. Additionally, 1 keV and 10~ cm" 
provide convenient fiducial cluster temperature and particle density values. In 
addition I assume the ICM is pure hydrogen and that protons and electrons 
share the same temperature. 



2. How vi^ell does MHD apply to the ICM? 

MHD is widely applied to the dynamics of the ICM. Being a single-fluid, contin- 
uum model, MHD requires on timescales of dynamical interest, ~ id/ud, that 
the particle populations pass through a sequence of local equilibria; i.e., effec- 
tive particle interaction times, tc ^td-, and interaction lengths, Ic ^ ^d-, where 
Id and Ud are representative dynamical lengths and speeds. In addition MHD 
requires relevant charged particles to be effectively magnetized; that is, their 
gyroradii, = vj_mc/{eB) <C id and gyroperiods, tg = 27rrg/v± = 2-ir/u;c td- 
MHD usually neglects the displacement current in Ampere's Law, a valid ap- 
proximation if relevant wave speeds are nonrelativistic; i.e., oj <C kc. Standard 
MHD also assumes isotropic (scalar) pressure and transport properties for the 
fluid, global charge neutrality. It neglects electrostatic forces due to local charge 
fluctuations, allowing the fields and bulk fluid to vary simultaneously. There are 
variants to the standard MHD model that relax various of these assumptions, 
such as a single fluid, p r essure and transport isotr o py (e.g., Bra-g inksii 196^; 
iBovd Sanderson 20031 : IGoedbloed Poedts 200i : iKulsrud 200il 'l. although 



they are more complicated to apply. 

The appropriate interaction lengths are central to this discussion. Be- 
yond their role in validating the continuum dynamics model itself they also 
control basic fluid properties needed in the MHD model, including transport 
behaviors such as viscosity, u, thermal conductivity, k, and electrical resistiv- 
ity, -q. T he standard kinetic theor y expressions for these transport coefficients 
are (e.g.. iBovd fc Sanderson 20031 ) v = kBTtc^p/m-p, k = [bkB) /{2me)nkBTtc^e 



and T] = rrieC^ /{ne'^tc,e), where the second subscript, 'p' or 'e' refers to protons 
or electrons. When tc is the Coulomb collision time discussed below the ther- 
mal conductivity is known as the "Spitzer conductivity", while the analogous 
viscosity is called the "Braginskii viscosity" . 

Unlike laboratory fluids where strong binary collisions typically establish 
equilibria, the hot, rarefled and fully ionized conditions in an ICM lead to a 
dominance by the combined effects of many weak interactions. The number of 
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charged particles participating in Coulomb collisions with a proton or electron 
is determined by the number of particles, Nd, found inside the so-called Debye 
sphere of radius, Xd, given by 

Ad = [A:bT/(4W)]1/2 = [^,\/(3^2)]i/2 ^ 2.4 x 10^ {Tk,v/n -2)'^^ cm, (1) 
where 

vth = ^JskBT/m = 5.4 X lO^T^/J {rrip/m)^/^ cm/s (2) 
is the characteristic particle thermal speed and 

ujp = sj47rne'^/m = 1.3 x 10^ {mp/m)^/^n]!2 s"^ (3) 

is the plasma frequency for that species. In the numerical expressions, nip is the 
proton mass and the plasma density is expressed in units n = 10~^n_2 cm~^. For 

electrons, vth = 2.3 x W^T^^y cm/s and u!p = 5.6 x lO^nV^^ s"-*^. When Nn » 1 
the Coulomb scattering cross section is enhanced over that for strong binary 
collisions, ac ~ Ane'^ /{SksT)'^ , by a factor proportional to ln(3A''£)) = InA due 
to random, thermal fluctuations in charge density. 

In clusters No = (47r/3)neA|, ^ 5.4 x 10^^{TkeV /n^2)^/^ > 1, and InA ^ 
30—40. Using In A = 35, the effective proton collision time is then ( Spitzer 19621 ) 



tc,p ~ 4.1 X lO^rf/^/n.s yrs, (4) 

while the electron collision time, fc,e, is smaller by a factor [me/nipY^'^ . This 
gives both particle species a collisional mean free path , £c = vthtc-, given by 

4,p,e ~ 22 Tlll/n.2 PC. (5) 

Temperature equilibrium between protons and electrons actually requires 
more stringent constraints, since the energy exchange during e-p collisions is pro- 
portional to nie/mp. Consequently the time for thermal equilibration between 
the two species is ~ {mp/meY^'^tc,p- In relatively denser and cooler ICMs these 
various Coulomb times and lengths should be short enough to support a fluid 
model on many scales of interest. However, in less dense environments outside 
cores and particularly in hotter ICMs Coulomb interactions can be uncomfort- 
ably slow in this context. With n = 10^^ cm^^ and T = IQkeV we have from 
equations S] and [5] tc,p ~ 10^ yrs and Ic ~ 10^ kpc, for instance. 

On the other hand, except for uniform, static and unmagnetized media. 
Coulomb collisions probably underestimate the effective interactions available. 
As a starter, velocity space plasma instabilities in a dynamical setting may re- 
distribute particle motions more rapidly than Coulomb collisions. For example, 
if there are density or velocity gradients on scales smaller than the Coulomb 
interaction lengths, particles should stream against these gradients. Then the 
so-called two-stream instability can come into play. The two-stream instability 
leads to particle bunching and associated coherent electrostatic fields on scales 
of the Debye length that can redistribute particle momenta. In the simplest 
case two equal, cold, like-charged beams interpenetrate with a relative speed Vh- 
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Fluctuations with w avenumber k ~ ujp/vi, are unstable with very fast growth 
rates F ~ a;„ fe.g- IStix 1991 ). If we associate Vf, with the eventual thermal 



speed, then kXij ~ 1, and we see that particle motions are redistributed very 
quickly on scales of the Debye length. More realistically only a fraction of the 
total particle population would be involved in streaming, so growth rates would 
be reduced proportionately. Still, the coherent nature of the induced charge 
fluctuations can substantially boost the effective scattering rate compared to 
incoherent Coulomb scattering. 

Magnetic fields obviously influence this picture substantially, as well, since 
proton and, especially, electron gyroradii, r^, and gyroperiods, tg = 2iT/ujg, in 
the ICM should be small; in particular, rg/£c <C 1 or equivalently, tcUJg S> 1. For 

thermal protons and electrons, respectively, Vg^p = 4.6 x l^fT^j^y/B^c cm, r^^e = 
1.0 X W^T^I^y/ Bfj^G cm, and with LOg = eB/{mc), ujg,p = O.OIB^g sec"-*^, tOg^e = 
ISBfj^G sec~^. Accordingly, 

tc,pi^g,p = {me/mp)'^/'^tc,ei^g^e ~ 1-2 X W^^T^^y B^cnZ], (6) 

which will generally be large in the ICM. This will guarantee, as well, that the 
ICM is magnetized sufficiently to apply MHD, whenever ic ^ ^d- 

MHD assumes charge quasi-neutrality, so that electrostatic fields can be ne- 
glected. In effect, when there is global neutrality one requires inside fluctuations 
that \q\/{ne) <C 1, where q is the local charge density. Dimensional analysis com- 
bining the equation for charge continuity, dq/dt = —V • j, with Ampere's law 
and the previously mentioned nonrelativistic phase speed constraint leads to the 
condition u}/u>g^e iy^p,el^g,e)^ ~ 10 n-2/B'j^Q, where the condition is applied 
to electrons because of their greater mobility. This limits us to ICM fluctua- 
tions roughly slower than a gyroperiod, so that free electrons can redistribute 
themselves to short out local electrostatic fields. 

In a quasiuniform field, particles spiral along the field with orbital radius, 
Tg, traversing a longitudinal distance, ic ~ v\\tc before scattering substantially 
changes their pitch angles. Scattering also introduces a transverse gyrocenter 
shift ~ Tg, so in the limit rg/ic ^ 1, the relative transverse and parallel particle 
diffusion rates would be ~ {rg/icY. 

Thermal energy and also momentum transport across fields are controlled 
by cross-field diffusion, which reduces by similar factors, ~ (r^/^c)^, the ther- 
mal conductivity transverse to a unifor m field and also the viscosity i n response 
to velocity gradients transverse to B ( Braginksii 19651 : ISpitzer 19621 ). The co- 



efficients parallel to B are also modified by dimensionalit y influences, as is the 
electrical resistivity tensor, by factors of order unity (e.g.. iKulsrud 2003 ). 

However, the ICM is thought to be turbulent w ith turbulent motions con- 
tributing perhaps ~ 10% of the total pressure fe.g.. ISchuecker et al. 20"oi ). In 
that case the magnetic field might not be at all uniform. Spatial diffusion then 
could be limited by the scale for bending of field lines or wandering of individual 
field lines. There has been much disc ussion of this topic recently, especially with 
regar d to thermal conduct ivity (e.g.. iNaravan &: Medvedev 200 ll ) and viscosity 



(e.g., Revnolds et al. 20051 ). Lazarian has recently provided a nice outline of the 



issues (jLazarian 20071 ). 



MHD Interactions 



5 



For our discussion the relevant points would be these. The character of the 
turbulence depends on whether U(i[t)/vA is greater or lesser than unity, where 

VA = S/(47rnmp)i/2 22B^GnZl^'^ km/s (7) 

is the Alfven velocity and Ud{i) is the turbulent velocity on scales, i. When 
Ud{t}/vA is large the turbulence is isotropic and Kolmogorov, so that Ud[tj oc 
£i/3_ "When U(i{(.)/vA is small, on the other hand, field line tension is sufficient 
to resist bending, so that turbulent motions and magnetic field structures be- 
come decidedly anisotropic. On the turbulence injection scale, L, in clusters 
the condition Ud{L)/vA ^ 1 probably applies under most circumstances. The 
transition to anisotropic, MHD turbulence then occurs below a scale 

iA/Lr^lQ-^Bla/^ulimnlh. (8) 

where ml,ioo is the injection-scale turbulent velocity in units of 100 km/s. Par- 
ticles should free stream no farther than the lesser of Ic and Ia- On scales larger 
than Ia the particle diffusion (and the associated transfer coefficients) should 
be isotropic with 

K ~ {l/'3,)min{l, ^A/ DuSpitzer < {l/3')KSpitzer (9) 

and 

ly ~ {l/5)min{l,iA/Ql^Braginskii < {'^ / 5)'^Braginskii (10) 

(e.g., iGruzinov 200l : |G ruzmov 20061 1. assuming Coulomb scattering. Prom the 



relations already given we can estimate that 

-5/2 -3/2 



- O.^B;GT,,;^nZ;'ul^,ooLkpc. (H) 

Por common ICM conditions IaI^c < 1 and sometimes IaI^c giving viscosi- 
ties and the thermal conductivities at least several times smaller and perhaps 
much smaller than Coulomb values. 

Pinally, standard MHD assumes an isotropic plasma pressure. To under- 
stand this constraint consider particles in a smoothly varying magnetic field. 
Without collisions the particles move adiabatically, conserving v\/B and, under 
incompressible conditions, v\\B. Again assuming incompressible changes and 
defining a pressure anisotropy, A = (Pj_ — P||)/(2P||), with P = {l /2>)nv^^ it is 
easy to show starting from ^ = that dA/dt ^ {3/2)dln{B)/dt (jHall 1980l ). 
This growth in pressure anisotropy is limited to the collision time. So, if we 
choose td for changes in the magnetic field, we have A ~ {3/2)tc/td- Thus, 
where collisions are fast on dynamical times the pressure anisotropics, A, should 
be small, and a scalar pressure representation adequate. Pirehose and mirror 
instabilities that develop easily on scales down to a few gyror adii from pressure 
anisotropics when magnetic fields are weak (e.g., iHall 1980l ) also will tend to 
reduce anisotropics, while enhancing magnetic field irregularities. 

We can summarize this section by saying that if we must depend entirely 
on Coulomb scattering, then interaction lengths may be sufficient in cooler and 
denser ICMs to support the MHD model, but not so in hotter and more rarefied 
ICMs. On the other hand plasma processes and even weak magnetic fields may 
substantially change that situation, allowing MHD to be a meaningful model at 
least beyond kiloparsec scales in most ICM environments of interest. 
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3. Ideal MHD? 

Accepting some caveats from the previous section, we can proceed and express 
the standard dissipative MHD equations in terms of an isotropic gas pressure, 
P, and scalar viscosity, v, thermal conductivity, k, and resistivity, r] in electro- 
magnetic units as 

§f + V • ipu) = 0, (12) 

+ li • Vu + - VP - - — (V X B)x B = i^V'^u, (13) 
ot p Attp 

^ + ^ . VP + 7PV • n = (7 - 1) [V • (kVT) + ^(V x B)% (14) 

^-Vx(.xP) = ^V^P, (15) 

V • P = 0. (16) 

The equations represent conservation of mass (eq. [T2]) . momentum (eq. 
[T3]) . energy (eq. and magnetic flux (eq. [T6]) plus Faraday's induction law 
combined with Ohm's law (eq. [T5|) . It is most common to apply the ideal, or 
nondissipative version of these equations, in which the viscosity, resistivity and 
thermal conductivity are neglected. How appropriate is ideal MHD in the ICM? 

Electrical resistivity is the easiest to evaluate. A comparison of the inductive 
and resistive terms in eq. [15] leads to the condition that resistive dissipation and 
field diffusion can be neglected when the magnetic Reynolds number, Rm = 
Ud(^d/v ^ 1- Assuming Coulomb collisions, 

Rm ~ 2.3 X 10^\d,iooid,kpcT^iy, (17) 

so that, as expected, the magnetic field is nicely frozen-in on most scales of 
interest. Ohmic heating can be neglected in eq. O when RmI^ ^ 1, where 
P = P/Pb, with Pb = PV(8vr). We have (5 ^ 400n_2rfcey/p2^ > 1, so that 
Ohmic heating is generally unimportant. 

Thermal conduction in eq. [TH can be neglected with respect to convective 
heat transport when the Peclet number, Pg = Udf-dP/ii^T) ^ 1- With Spitzer 
conductivity we have n 2.6 x lO^^T^^y erg/(cm s K), so that 

Pe ~ 3.3 X 10-\d,ioo^d,fepcn-2Pfcey'. (18) 

Fig. 1 illustrates Pr from this expression. It is clear, as many have noted before 
(e.g., Bertschinger fc Meiksin 198^ : iNaravan k, Medvedev 200il : iFabian et al. 2005 



Pope et al. 20051 ) that Spitzer conductivitv would not be ignorable in the ICM 



However, as discussed above, especially if the magnetic field is tangled on scales 
small compared to the Coulomb scattering length the effective conduction could 
be substantially smaller than Spitzer. How small remains an open question. 
Cold fronts, which are contact discontinuities sometimes seen to exhibit sharp 
temperature changes, even on scales less than a Coulomb scattering length 
( Vikhlinin et al. 200ll ). demand a much smaller conductivity, at least locally. 
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The importance of viscosity is mea- 
sured by the relative size of the inertial 
to the viscous terms in eq. [13l that is, 
by the Reynolds number, Re = Udid/i^- 
When Re ^ 1 inviscid dynamics is 
a good approximation. The Reynolds 
number must also be large for turbu- 
lence to be established for scales, id and 
Ud- As for thermal conduction, it has 
been pointed out that Coulomb (Bra- 
ginskii) viscosity can sometimes be too 
large in the ICM to assum e these limits 
(;e.g.. lRevnolds et al. 20051 '). Here 

Re ^ lSud,imU,kpcn-2Tj^eV^ ■ (19) 

Re and Pe are related by the Prandtl 
number, Pr = Pel Re ~ 0.08(me/mp)V2^ 
reflecting the ion diffusion and electron 
diffusion origins for viscosity and ther- 
mal conduction. 

Fig 1. illustrates the Braginskii ex- 
pression for Re- As before, in cooler, denser ICM environments this estimate 
often gives i?e ^ Ij whereas in hotter rarefied ICMs that is not the case. One 
would have to reduce the viscosity two orders of magnitude below the Bragin- 
skii value to establish < 100 km/s turbulence on 10 kpc scales if T^eV = 5 and 
n_2 = 0.1, for instance. However, if turbulence is somehow established down to 
scales below ic{C oulomh) , perhaps initially aided by plasma processes, the effec- 
tive particle streaming lengths may become small compared to allowing the 
large Re needed for turbulent conditions to be maintained. Given observations 
that suggest turbulence in ICMs on scales of a few kpc and recognizing that ra- 
dio halos may depend on tu rbulent cosmic ray reacceleration from small eddies 
( Brunetti &: Lazarian 200?! ) it is perhaps reasonable to apply inviscid MHD on 
those scales. We may hope that ultimate clarification of this issue will come from 
definitive measures of turbulence in a wide variety of ICM settings, especially 
on the smallest detectable scales. 

4. MHD influences on AGN/ICM interactions 

I assume thus, with some reservations, quasi-ideal versions of eqs. [T2]|151 Our 
primary interest, in any case, is the dynamical influence of the field, as expressed 
by the Lorentz force term in eq. [T3j The standard criterion for evaluating its 
importance is to compare the magnetic pressure portion of the Lorentz force to 
the gas pressure gradient, so that magnetic effects are assumed small whenever 
j3 » 1. In the ICM we have /? « AQ{)n-2TkeV / B'j^Q. Using this criterion alone 
we would conclude that Maxwell stresses in the ICM were unimportant. That 
is probably a false impression, however. 

While a meaningful rule of thumb, the simple (3 metric overlooks important 
details. A particularly simple one is that the magnetic and gas pressure forces 




0.0001 ' ' ' — ' — ' ' ' ^ ^ 

1 10 
T 

' keV 

Figure 1. Reynolds number, 
Re, (upper, black lines) and 
Peclet number, Pg, (lower, red 
lines) assuming Braginskii vis- 
cosity and Spitzer conductivity. 
Densities are n — 10~^ cm~^ 
(solid) and n — 10^'^ cm^^ 
(dashed) . 
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depend on is and ip, the spatial variation scales for B and P. These scales 
need not be the same, especially in a complex, mostly incompressible, turbulent 
setting. Then, when /?>>!, the magnetic field, being tangled and intermittent 
in magnitude, may exert a force that is much stronger than one would anticipate 
from a simple evaluation of (3. A better met ric would be I3{Ib/^p) ^ 1- This 
point has been made by (lO'Neill et al. 20071 ^ in a 3D MHD simulation study 
of AGN jets in an ICM containing a tangled magnetic field. Even for a mean 
f3 ~ 100, they found that the ICM magnetic field had a very significant influence 
on the subsonic expansion of the jet/ICM contact discontinuity, because local 
values of (3 varied by factors of several on scales that were much smaller than the 
relatively smooth ICM gas pressure. Consequently magnetic pressure gradients 
were, in fact, sometimes comparable to gas pressure gradients. 

A second point is that the magnetic tension force, which resists bending 
of field lines, is often more important than the pressure force, especially in 
subsonic flows characteristic of ICMs and many of their interfaces with AGN 
outflow residues. An example setting where this applies, mentioned in §2, is 
MHD turbulence, where magnetic tension distorts small scale eddies. The role 
of magnetic tension can be estimated by comparing the Maxwell stress in eq. 
[T3] to the inertial, shear stress. From this we establish that magnetic tension 
may be important when va > U(i{^b / ^d)- If ~ ^d, as is the case parallel to 
the field in MHD turbulence, we recover the va ^ u criterion from §2. Other 
examples will be discussed below. 

Another point to keep in mind, of course, is the possibility of magnetic 
field amplification, described by the induction eq. [13l Outside reconnection 
regions where dissipation is clearly important, the ideal MHD evolution of the 
field expressed in eq. [T3] consists of field compression and field line stretching. 
In complex flows field line stretching is usually a bigger effect, since it is limited 
only by the eventual backreaction of magnetic tension and local reconnection 
when field lines "cross". For incompressible flow the ideal version of eq. [T3l 
reduces to \B\ = |i?o|^/^0; where £ measures the length of a field line segment. 
Thus, the field intensity is proportional to the stretching of the line. This will 
vary stochastically in a turbulent fiow, leading to spatially intermittent field 
intensities. The constraint imposed by magnetic tension, va < on the other 
hand, assures that amplification of the magnetic field will be limited to magnetic 
pressures, Pb = {l/2)v\p, less than the turbulent, kinetic pressure, Pk ~ u'^p- 

One very important ICM MHD concern is the evolution of hydrodynami- 
cal instabilities, particularly the Kelvin-Helmholtz (K-H) and Rayleigh- Taylor 
(R-T) instabilities. These are both essentially incompressible instabilities that 
occur along boundaries such as the contact discontinuities separating AGN jet 
cocoons from the ICM or relic radio bubbles from the ICM. Related, R-T-like in- 
stabilities can develop anywhere equilibrium pressure and density gradients take 
opposite signs, or equivalently where pressure and entropy gradients have the 
same sign. It turns out that surprisingly weak magnetic fields can significantly 
modify these instabilities, sometimes quenching them and sometimes control- 
ling their nonlinear evolution even when magnetic stresses do not prevent the 
instabilities from developing. 

The K-H instability develops hydrodynamically along a shear layer for wave- 
lengths exceeding the shear layer thickness, provided the shear velocity is not 
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highly supersonic. Tension from a magnetic field component parahel to the ve- 
locity in the shear layer and aligned with the wavevector of a perturbation will 
resist growth of the instability, since growing oscillations lengthen the bound- 
ary. For a velocity discontinuity, Am, separating fluids of e qual density the linear 
growth rate of a small perturbation with wavenumber k is ( Chandrasehkar 196ll ) 



r = (l/2)A;AM^(l-4t;2 ||/An2). (20) 
The perturbation is stable (F^ < 0) for iiyi || > An/2. Linear stabilization of 

1/2 

the K-H instability in the ICM would roughly require B^g > ^n^^ ioo't-_-2 ' ^ 
condition that is likely to be met along interface s of low density and modest 
velocity contrast, as pointed out, for example, by ( De Young 20031 ). 



It is important to realize, however, even when it is initially much too weak 
to inhibit the linear K-H instability, that a magnetic field can play a very impor- 
tant role in the nonlinear K-H instability. This influence results from stretch- 
ing (amplificatio n) of the initial field within the unstable fiow. For example, 
( Rvu et al. 2000l ) carried out 3D MHD simulations of K-H unstable flows, and 



found that even for initial conditions with || ~ 0.02An, the unstable fiow 
evolved into a relaxed, laminar form with aligned, self-organized magnetic and 
velocity fields, contrary to the turbulent velocity field that developed in the ab- 
sence of a magnetic field. The key insight is that the 3D hydrodynamical flow 
initially forms into line vortices along the unstable slip surface, stretching the 
magnetic field across each vortex. Subsequently, the vortex is unstable in 3D to 
so-called rib vortices transverse to the main vortex. This greatly adds to the net 
stretching of the field, making it possible locally and temporarily that va > An 
and, thus for an initially weak field to alter evolution of the fiow. 

The linear R-T instability is infiuenced by a magnetic field parallel to a 
density discontinuity in a way similar to the K-H instability; magnetic tension 
resists the stretching of the boundary. Defining the density jump across the 
discontinuity as r = P2/ pi, where region '2' is on top with respect to gravity , g, 
and letting v\ = i?^/[47r(pi + P2)], the growth rate is ( Chandrasehkar 196ll ) 



[l/{r + l)]J[gk{r-l)-2vl,,k^]. (21 



In this case the influence of the magnetic tension depends on wavelength, so that 
the instability is suppressed when k > ko = {l/2)g{r — l)/[v^ ||(r + 1)], while 

there is a maximum linear growth rate at wavenumber, kmax = (l/2)A;o, given 
by 

Tmax = (l/2V2)<7(r- l)/[^^,||(r + l)]. (22) 

As a mnemonic, imagine for a large density contrast, r ^> 1, that the free fall 
over a time ~ 1/Tmax reaches a velocity ~ ^^,11) and a length ~ 1/kmax- It turns 
out even a vertical magnetic field reduces linear growth of the R-T instability, 
although it does not totally quench it. The asymptotic short wavelength growth 
rate is of order q/va, compared to the hydrodynamical rate ~ \/^- These two 
rates are similar for k ~ mirroring the horizontal field behavior. As a 

rule of thumb, then, we can expect magnetic inhibition of the R-T instability 
whenever the Alfven speed is greater than or comparable to yfgl- That would 
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correspond roughly to a magnetic field B^g > 7((n_2 g-Tid.kpcY^"^ ■, where g--j is 
the gravitational acceleration in units of 10~^ cm/s~^. Alternatively one could 
view this relation to indicate, given ICM magnetic fields of order a microGauss, 
that R-T instabilities on scales smaller than a few kpc ar e likely to be inhibited . 
As I noted earlier and others have also emphasized (e.g., iRevnolds et al. 20051 ). 
viscous effects on this scale could also be important, depending on how much 
plasma and magnetic fields reduce the effective free-streaming length for protons. 
Which influence, magnetic tension or viscosi ty, is dynamically more important 
on this scale remains to be established (e.g., iKaiser et al. 2003 ). 

Ideal MHD influence s in the nonline ar R-T instability are a bit compli- 
cated, as pointed out by ( Jun et al. 19951 ) in a 2D and 3D simulation study of 
the nonlinear MHD R-T instability. The nonlinear hydrodynamic R-T instabil- 
ity is characterized by the formation of dense fingers that 'drip' downward and 
light "bubbles" that rise. These eventually lead to turbulent conditions. Jun 
et al. noted that the vertical component of B aligned with gravity eventually 
dominates the evolution of these structures, especially when the initial field is 
only moderately weaker than needed to suppress the linear instability. Then, 
in fact, since the vertical magnetic field aligns with the edges of R-T fingers, 
it inhibits development of a secondary R-T instability that otherwise disrupts 
those fingers. Thus, it can actually enhance their early development. Ultimately, 
however, even a relatively weak magnetic field caught up in the nonlinear de- 
velopment of the R-T instability can be amplified through stretching to restrict 
growth of turbulent behaviors. 

There are a limited number of MHD simulation studies of AGN and/or 
radio relic, bubble interactions with ICMs. They generally support the points 
made here, including the realization that fields at least several times weaker 
than those expected by the usual metri cs can play essential dynamical roles. As 
examples, ( Robinson et al. 20041 ) and ( Jones &: De Young 2005 ) considered 2D 
MHD buoyant bubbles in model ICMs with large scale magnetic fields, point- 
ing out that even when (3^1 field stretching could st abilize boundaries that 
otherwise were disrupted by R-T and K-H instabilities. ( De Young et al. 20071 ) 
obtained similar results for 3D bubbles, but also noted that since field line ten- 
sion acts only in the plane containing the field and its curvature, the bubbles 
remained subject to disr uption along lines orthogonal to the field-gravity plane. 
( Ruszkowski et al. 20071 ) carried out 3D MHD bubble simulations in which the 
ICM magnetic field was tangled. They noted that when the outer tangling scale 
was smaller than the size of the bubble the magnetic fiel d was no longer effectiv e 
in stabilizing the bubbles. This result, confirmed by ( De Young et al. 20071 ). 
reflects the fact that disruption comes from instabilities on the scale of the 
bubble. It is obvious that magnetic tension, the primary MHD stabilizing in- 
fluence, does not operate on scales greater than the coherence length of the 
field. On those scales we should expect the dynamics to be largely hydrody- 
namic, except for m agnetic pressure gradient effects. I already mentioned that 
( O'Neill et al. 20071 ) pointed out how variations in magnetic pressure on scales 
smaller than gas pressure variations can produce important dynamical conse- 
quences even when (3^1. 
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5. Conclusion 

As a number of previous authors have discussed, Coulomb mean free paths in the 
ICM can sometimes be uncomfortably large when one wants to model the ICM 
as an ideal fluid. Yet observed features, such as shocks, sharp contact discontinu- 
ities and turbulence strongly suggest continuum, fluid behaviors. It seems likely 
that a combination of plasma instabilities and magnetic field influences reduce 
the effective particle free-streaming lengths enough to allow relatively inviscid, 
fluid-like behaviors at least beyond kiloparsec scales. Whether these effects can 
also reduce diffusion rates of electrons sufficiently to effectively quench thermal 
conduction is less certain. In any case electrical conductivity is very likely to be 
sufficient to apply a frozen- in magnetic field model on these scales. 

Even though magnetic pressures in the ICM are likely to be much smaller 
than gas pressures, magnetic fields there can still play a significant dynamical 
role. This is especially because the local Alfven velocity can be comparable to 
or exceed local flow speeds, since ICM flows are usually subsonic. In that case 
magnetic tension forces can compete with dynamical stresses. This probably 
influences ICM turbulence on small, say kiloparsec, scales and helps stabilize 
flows that are otherwise hydro dynamically unstable. 
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